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ABSTRACT The recent BeppoSAX observations of the BL Lac objects Mkn 501 and 1ES 
2344+514 have shown that the synchrotron spectrum of these objects peaks, in a V-vF v repre- 
sentation, at energies at or above 100 keV. One can wonder if these are the most extreme examples 
of hard synchrotron blazars, or if they are the first cases of a more numerous class of sources. Here 
I propose the existence of a class of even more extreme BL Lac objects, whose synchrotron spec- 
trum peaks at or above 1 MeV. Based on the observational trend found between the location of 
the synchrotron peak and the bolometric power of BL Lac objects, it is argued that the proposed 
extreme sources could have escaped detection (in any band) so far, or could have been classified 
as galaxies, and their "BL Lac-ness" could be revealed by INTEGRAL. 

KEYWORDS: BL Lacertae objects; synchrotron emission; inverse Compton emission; radio jets; 
X-rays and gamma-rays: spectra 

1. INTRODUCTION 

The blazar class of sources is formed by BL Lacertae objects (BL Lacs) and flat 
spectrum radio quasars (FSRQ), either with high (HPQ) or low (LPQ) polarization 
in the optical band. The differences among subclasses of blazars reflect the presence 
or absence of emission lines, the degree of optical polarization, the band (i.e. radio 
or X-rays) where they were discovered. 

As an example, consider BL Lacs discovered through radio and X-ray surveys: 
they show different radio to X-ray spectra, but they share other properties such as 
the absence of strong emission lines, the rapid and large amplitude variability and 
the same average X-ray luminosity. This led Maraschi et al. (1986) and Ghisellini 
& Maraschi (1989) to propose that the spectral differences were due only to the 
different viewing angle under which we observe an accelerating, inhomogeneous jet. 

On the other hand, Giommi & Padovani (1994) noticed that the spectral energy 
distribution (SED) of radio and X-ray selected BL Lacs showed peaks (in a v-vF v 
representation) at different energies, and suggested that this difference was intrinsic, 
and not due to orientation effects. They therefore introduced the notation of HBL 
(high energy peak BL Lac) and LBL (low energy peak BL Lac), the former being 
sources preferentially selected through X-ray surveys, and the latter through radio 
surveys. A crucial help to understand blazars came from EGRET, onboard CGRO, 
and from the ground based Cherenkov telescopes, such as Whipple and HEGRA: 
we now know that most of the power emitted by blazars often lies in the 7-ray 
region of the spectrum. Their SED presents two peaks, the first at mm to X-ray 
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energies, and the second in the MeV to the TeV band. The energies of the two peaks 
correlate, in the sense that if the first peak is located in the mm band, the second 
is at MeV energies, while if the first peak is in the X-ray band, the second is at 
GeV-TeV energies. The idea of Padovani & Giommi (moving peak) can therefore 
be extended also to the second peak. 

Understanding what rules the SED of blazars is one of the main goal of the 
current blazar research. This is not an easy task, however, since there is still dis- 
cussion about the nature of the radiation we see: while the emission from the radio 
to UV (i.e. the first peak) should be due to the incoherent synchrotron process, the 
emission at higher energies (second peak) could be pure Synchrotron Self Compton 
(SSC: Bloom & Marscher 1996; Maraschi, Ghisellini & Celotti 1992) or a mixture 
of SSC plus a contribution by inverse Compton scattering off photons produced ex- 
ternally to the jet (EC: Sikora, Begelman & Rees 1994; Dermer & Schlickeiser 1993; 
Blandford & Levinson 1995, Ghisellini & Madau 1996), or another more energetic 
synchrotron component (as in the 'proton blazar' model by Mannheim 1993). In 
this framework, objects presenting extreme properties are the most interesting ones, 
since they best constrain our models. In this respect the discovery made by Bep- 
poSAX of an extraordinary hard X-ray state of the BL Lac objects Mkn 501 and 
1ES 2344+514 is extremely important, and stimulates new and interesting ideas on 
the physics of relativistic jets. 

In this paper I will review some recent work on the SED of blazars, both from 
the phenomenological and from the theoretical point of view. Then I will discuss 
the possibility that an entire new class of BL Lacs can exist, characterized by a 
synchrotron component (i.e. the first peak of the SED) peaking at MeV energies. 

2. THE BLAZAR SEQUENCE 

Fossati et al. (1998), by considering flux limited samples of BL Lac objects and 
flat spectrum quasars, have noted an interesting relation between the SED of blazar 
and their bolomctric observed luminosity. This is shown in Fig. 1. In less powerful 
objects the synchrotron peak is at higher energies, reaching the soft-medium X- 
ray range. Analogously, the high energy peak, believed to be due to the inverse 
Compton process, shifts to lower energies as the total power increases. In addition, 
the ratio between the 7-ray power and the optical-UV luminosity increases as the 
total power increases. 

Ghisellini ct al. (1998), after modelling all blazars detected by EGRET with 
an homogeneous EC model, interpreted these results on the basis of the underlying 
correlation between the Lorentz factor of the electrons emitting at the peak, 7 pea k) 
and the amount of energy density (both magnetic and radiative) present in the 
emitting region (see Fig. 2). This in turn correlates with the observed (beamed) 
luminosity, and the ratio between the power of the Compton to the synchrotron 
components. Blazars form a sequence: low luminosity lineless BL Lacs (HBL) 
have large values of 7 poa k, synchrotron peak energy in the EUV-soft X-rays and a 
roughly equally powerful Compton component peaking in the GeV-TeV band. LBL 
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FIGURE 1. The average spectra of blazars, according to their bolometric luminosity, assumed 
to be traced by the radio luminosity. Blazars belonging to complete samples have been devided 
in 5 luminosity bins, irrespective of their classification. Notice how the average SED changes 
as the overall power changes. Dashed lines corresponds to an analytical, phenomenological fit 
to the data. From Fossati et al. 1998. 
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FIGURE 2. Correlations found when modelling the SED of EGRET blazars with an ho- 
mogeneous EC model. The left panel shows the correlation between "fp ea k an d the total 
(magnetic plus radiative, including the contribution from external photons) energy density. 
Different symbols indicate different classes of blazars, as labelled. The right panel indicate the 
correlation between the compactness [t EE L(Tt / (Rm c C^) of the external photons and the 
one corresponding to the injected power, as measured in the comoving frame of the emitting 
region. From Ghisellini et al. 1998. 
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are characterized by a greater overall luminosity, a smaller 7 pea k and peak energies 
in the optical and GeV band. More powerful sources, such as HPQ and LPQ, have 
the smallest values of 7 pea k, peak energies in the mm-far IR and MeV band, a 
dominating Compton component in which photons produced externally to the jet 
are more important (as seed photons to be Compton scattered at high energies) 
than the locally produced synchrotron photons (see Fig. 2). 

2.1 The case of Mkn 501 and 1ES 2344+514 

In April 1997, BeppoSAX observed Mkn 501 (Pian ct al. 1998), one of the closest 
BL Lac objects (z = 0.034), and the second extragalactic source detected in the TeV 
energy range by Cherenkov telescopes (Quinn et al. 1996; Bradbury et al. 1997). 
Results of the X-ray observations were extraordinary (Pian et al. 1998, see also Pian 
et al., these proceedings), since BeppoSAX detected the source in an extremely high 
and hard synchrotron state, simultaneously with a TeV flare. The vF v synchrotron 
spectrum of this source was observed to peak at 100 keV or beyond, while, at 
~0.5 TeV, Mkn 501 was a factor 4-6 brighter than the Crab (Catanese et al. 1997, 
Aharonian et al. 1997). A similar hard X-ray spectrum was observed by BeppoSAX 
also for another TeV source, 1ES 2344+514 (Giommi, Padovani & Perlman, 1998). 
Extreme objects like Mkn 501 and 1ES 2344+514 could be sources where the particle 
acceleration mechanism operates at its maximum efficiency, succeding to accelerate 
electrons up to 10 TeV or more. 

Mkn 501 and 1ES 2344+514 are probably in the extremely hard state observed 
by BeppoSAX only when flaring. However, there could be other sources that are 
extremely hard even in quiescence. The above discussed trends suggest that these 
sources should be at the lower luminosity end of the BL Lac luminosity function. 
It is interesting then to discuss what can limit the relevant electron energy 7 pca k, 
which in turn determines where the synchrotron spectrum peaks. 

3. LIMITS TO THE ELECTRON ENERGIES 
3.1 Shocks 

Guilbert, Fabian & Rees (1983) derived a useful limit on the maximum synchrotron 
frequency that can be produced by shock-accelerated electrons. In rclativistic 
shocks the Lorentz factor fractional change of the electrons for every passage through 
the shock can be of order unity (A7/7 ~ 1), with the acceleration timescale ap- 
proximately equal to the gyroperiod (oc 7_B _1 , where B is the magnetic field). 
The maximum energy 7 max is attained when the synchrotron cooling timescale 
oc 7~ 2 £? -2 equals the acceleration timescale. Then 7 max oc B -1 / 2 , resulting in a 
^-independent maximum synchrotron frequency of 70 MeV. Additional Compton 
losses would of course lower this value, but not severely, at least in a pure SSC 
model, since they are inhibited by the decline with energy of the Klein-Nishina 
scattering cross section. 
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With a Doppler factor (5^10 the observed maximum frequency could be ~ 700 
MeV. Values of this order lie within the EGRET capability. Therefore it is unlikely 
that BL Lacs have synchrotron spectra reaching energies greater than ^100 MeV. 
There must then be a more severe limit to the maximum observable synchrotron 
frequency. 



3.2 Global energetics 

Another limit to the acceleration of particles may come from the total available 
power. This is likely in the form of bulk kinetic motion of the emitting plasma, 
Lk, and in Poynting flux, Lb- They are given by (see e.g. Celotti & Fabian 1993, 
Ghisellini & Celotti 1998): 

L k = TtR 2 T 2 f3cri m c c 2 (< 7 > +m p /m ); L B = Ir 2 T 2 [3cB 2 , (1) 

8 

where R is the cross sectional radius of the jet, n' = Tn is the comoving particle 
density of average energy < 7 > m c c 2 , and m p , m c are the proton and electron 
rest masses, respectively. An electron proton plasma is assumed. The synchrotron 
intrinsic power is 

L' s = Volume J n'(7)7 s m c c 2 d-y = R 3 <r T cn'B 2 < 7 2 > (2) 

where % is the synchrotron cooling rate, and ^'(7) oc 7 _p between 7 m i n and 7 pea k- 
For a viewing angle ~ 1/T, the luminosity calculated assuming isotropy is related to 
L' s by i S! obs = r 4 is- The intrinsic power emitted over the entire solid angle equals 
T 2 ig. We can then relate the synchrotron power T 2 L' S to Lk (which is proportional 
to n') and Lb (which is proportional to the magnetic energy density Ub), obtaining 

Y 2 L' = — aTLkLB < 7 2 > 
s 9tt Rm c c 3 T 2 < 7 > +m p /m e 

Requiring T 2 L' S < Lk implies: 

< 7 2 > 9tt Rm c c 3 T 2 

■ < (4) 

< 7 > +m p /m c 16 <r T L B 

At high energies, the synchrotron process is efficient and fast, and a quasi-steady 
emission requires continuos acceleration of particles, at a rate that balances the 
radiative (synchrotron) losses. As discussed by Ghisellini & Celotti (1998), the 
most efficient synchrotron emitting jet corresponds to equipartition between bulk 
kinetic and Poynting powers. The observed synchrotron power is then maximized 
for Lk ~ Lb- Setting r 2 L^ at its maximum possible value (r 2 L^ ~ Lk ~ Lb), 
assuming n(-f) oc 7~ 2 between 7 m ; n < m p /m c and 7 pca k we have a limit for 7 pea k: 

9ir Rm u c 3 1 „ „ « Rm 
7 P cak < - a P — 7T = 1.2 x 10 6 fr- (5) 

a T 7min^ s .42 



5 



Here the notation Q = 10 X Q X is used, with cgs units. The observed synchrotron 
peak frequency ^ pca k = 3.7 x 10 6 7 2 cak -Br Hz corresponds to 

*Vak = 0-36 2 R lf\ 3/2 MeV (6) 

If only a fraction a of the bulk motion power is transformed into T 2 L' S , then the 
above estimate must be multiplied by a 7 / 2 . From Eq.(6) we have that the 1997 
flare of Mkn 501, with L' s ~ 2 x 10 42 /r 4 erg s" 1 and ^ pcak -100 keV, was quite 
close to transform all the available bulk motion power into synchrotron emission. 
For r = 15 and R = 10 16 cm, its emitted synchrotron power is ~40% of L k - 

The above arguments suggest that —1 MeV can be considered a limit for the 
observed maximum frequency of BL Lacs. If a is relatively constant, we also have 
an inverse correlation between f pea k and L' al in the sense that the most extreme 
objects in terms of ^ p0 ak are the least powerful ones. 



4. PREDICTED SPECTRA 

Figs 3 and 4 illustrate the possible spectrum of a MeV BL Lac, assuming as a 
working hypothesis an homogeneous one-zone SSC model. Rclativistic electrons 
are assumed to be continuously injected throughout a spherical source of size R 
embedded in a tangled magnetic field B. The injected particles have a power law 
energy distribution oc between 7 m ; n and 7 max . The steady emitting particle 
distribution is found by solving a continuity equation, i.e. by balancing injection 
and radiative cooling. The Klein-Nishina decline of the scattering cross section and 
the possible photon-photon pair production are accounted for. Details on the 
model can be found in Ghiscllini (1989) and Ghiscllini ct al. (1998). 

Also shown in Fig. 4 is the spectrum of Mkn 501 of the flaring state of April 
1997, together with other, non simultaneous, data. The computed SSC models 
obey the 7 poa k oc L^ 1 relation of Equation (6). The first of these models (top 
curve) assumes 7 max = 3 x 10 6 , L s ^ ohs = 2 x 10 46 erg s^ 1 R = 10 16 cm, 6 = 15 
and B = 0.9 Gauss. The intrinsic injected power is L( n j = 5.5 x 10 41 erg s _1 . The 
second model (second line from the top) is the same, but the assumed redshift is 
z = 0.1, and this redshift is assumed also for the remaining models. The size and 
the beaming factor are the same for all models. The magnetic field scales as L 1 / 2 , 
and each intrinsic (as well as observed) luminosity differs for a factor 3, with 7 max 
varying accordingly. All models assume I bs(v bs) = 5 3 I'(v r ) for the transformation 
between the observed and the intrinsic intensity. 

It can be seen that the plotted SED are at or below the current limits of large 
area X-ray survey, below the approximate sensitivity limit of EGRET, and the 
expected radio flux is at level of 1 mJy or less. But in the MeV region these sources 
are bright, even if currently below the COMPTEL capabilities. Note also that the 
Compton flux is severely inhibited by Klein-Nishina effects. However, the Compton 
flux is calculated assuming scattering between electrons and photons produced by 
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FIGURE 3. Some synchrotron self Compton spectra obtained by continuously injecting elec- 
trons up to different 7 ma x, as labelled, corresponding to different total luminosities. From 
Ghisellini et al., in preparation. 
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FIGURE 4. Same spectra as of Fig. 3, here shown in the V—vF v representation, assuming 
a redshift Z — 0.1 for all models but the top one, which is the same model as the second 
one, but for Z — 0.034 (redshift of Mkn 501). The dashed horizontal lines marks the level 
of 10 mjy at 5 GHz, approximately the level of the faintest X-ray selected BL Lacs, and of 
5 X 10 — 12 erg s — 1 between 0.3 and 3.5 keV, approximately the limit flux of the Einstein 
SLEW survey. The SED of Mkn 501 is shown for comparison. From Ghisellini et al., in 
preparation. 
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the same electron populations; any additional mm far infrared photons field can 
increase the emitted high energy flux. 

Note, from Fig. 4, that the 0.1-10 GeV flux is below the current EGRET 
sensitivity level, and that, in the optical band, the contribution of an underlying 
galaxy could hide the non-thermal continuum. In the X-ray band the flux could 
be strong enough to let these sources be present in moderately deep X-ray surveys, 
such as the Rosat RASS and the Einstein EMSS. However, note that the radio flux 
of these sources is at the mJy level, at or below the limit of the present large area 
radio sky surveys. This may be the reasons why we have not yet discovered them, 
i.e. some of these sources are too faint at all frequencies but the MeV band, and 
even if bright enough in the X-ray band to be included in current samples, they 
could have been classified as normal (radio-weak) elliptical galaxies. INTEGRAL 
could instead discover the brightest objects through their intense MeV emission. 
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